Are pinholes the cause of excess current in superconducting tunnel junctions? A 
study of Andreev current in highly resistive junctions 
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In highly resistive superconducting tunnel junctions, excess subgap current is usually observed 
and is often attributed to microscopic pinholes in the tunnel barrier. We have studied the subgap 
current in superconductor-insulator-superconductor (SIS) and superconductor-insulator-normal- 
metal (SIN) junctions. In Al/AlO^/Al junctions, we observed a decrease of 2 orders of magnitude 
in the current upon the transition from the SIS to the SIN regime, where it then matched theory. In 
Al/AlOrr/Cu junctions, we also observed generic features of coherent diffusive Andreev transport in 
a junction with a homogenous barrier. We use the quasiclassical Keldysh-Green function theory to 
quantify single- and two-particle tunneling and find good agreement over 2 orders of magnitude in 
transparency. We argue that our observations rule out pinholes as the origin of the excess current. 

PACS numbers: 74.45. +c, 74.40.Gh, 74.55. +v, 85.25.Am 



Superconducting tunnel junctions have become in- 
creasingly important devices in applications ranging from 
medical and astrophysical sensors to quantum comput- 
ing. A hallmark of tunnel junctions, key to these ap- 
plications, is their minimal dissipation. This dissipation 
is often parameterized by a subgap conductance G sg in 
parallel with an ideal tunnel clement. Even though G sg 
is relatively small at low temperatures, it is often ob- 
served to be orders of magnitude larger than what is pre- 
dicted by theory. This excess dissipation is emerging as 
a potential limitation in a host of new applications. For 
instance, this may be a source of energy relaxation in 
superconducting qubits [1] and in tunable resonators [2] • 
In single-electron turnstiles, this leakage may limit the 
ultimate accuracy of a future current standard [5]. 

To elucidate the problem, we note that the tunnel 
model of superconducting junctions predicts that G sg 
should decrease exponentially as a function of tempera- 
ture [3]. This is true for both superconductor-insulator- 
superconductor (SIS) and superconductor-insulator- 
normal-metal (SIN) junctions. In experiments with 
highly resistive junctions, however, an exponential de- 
pendence is observed for temperatures down to approxi- 
mately 10 % of the critical temperature of the supercon- 
ductor, after which a saturation is observed at a value 
a few orders of magnitude smaller than the normal-state 
tunnel conductance, Gm [HIS]. 

In spite of being universally observed, the origin of this 
temperature independent subgap current has remained a 
puzzle for decades. One of the accepted explanations is 
provided by multiparticle tunneling [7] and multiple An- 
dreev reflections [8J. This has been successfully used to 
explain the subgap current in low resistance junctions 
[5]. However, even from the first observation [TO], ex- 
periments on highly resistive junctions have revealed a 
drastic discrepancy between the measurements and the- 
oretical predictions. According to theory [?1[TTJ[T2], the 



subgap current in uniform SIS junctions at low temper- 
atures should have a series of current steps at voltages 
where different multiparticle processes are activated. In 
theory, the ratio of the current below and above each step 
is proportional to the junction transparency, V . There- 
fore, in junctions with T ~ 10 — 10 -6 one expects a 
current step near eV = 2 A with ratio 10~ 5 — 10~ 6 . How- 
ever, the experimental values typically saturate around 
10" 2 - 10" 3 [H [HI H3HIS1 - In SIN junctions, one simi- 
larly expects the ratio G S9 /Gn to be of order 1CT 5 — 10~ 6 . 
However, a similar saturation is typically observed |17l - 
|2"U] (though a smaller G sg was recently reported in [2"T]). 
This discrepancy has often been explained [7] by assum- 
ing microscopic defects in the tunnel barrier, commonly 
known as "pinholes," which cause a greatly enhanced lo- 
cal transparency (r > 10~ 3 ). In fact, a large subgap 
conductance has been considered an indicator of a poor 
quality tunnel barrier [3j l9l IT9] . 

In this Letter, we report on an extensive study of the 
differential subgap conductance G sg in tunnel junctions 
for a wide range of Gat. First, we have investigated 
Al/AlO^/Al junctions (Fig. [I]), similar to those used in 
qubit circuits. We investigated both the SIS and SIN 
transport regimes in these junctions. This was done by 
making one electrode thicker and applying an external 
magnetic field to suppress the superconductivity in this 
electrode. The two-particle current should not change 
significantly during this crossover according to theories 
for diffusive SIN [HNH] and SIS [H] junctions. However, 
in drastic contradiction to the theoretical expectation, 
we observed a decrease of up to 2 orders of magnitude in 
the subgap current in the SIN regime (Fig. [I]). G sg /Gjv 
then reaches values predicted by the theory of diffusive 
Andreev transport |12l 1221 123) , exhibiting a linear depen- 
dence on Gtv over more than 2 orders of magnitude. We 
ultimately achieve G S9 /Gn ~ 1CT 5 (Fig. [2]). 

To corroborate that we reached the fundamental limit 
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FIG. 1: (Color online) I-V characteristics of two typical low- 
transparency Al/AlOz/Al junctions (red and black) at T < 
100 mK. In the SIS regime (left, B = mT), a step is observed 
with the ratio of currents below and above eV = 2A being 
2 x 10~ 3 , consistent with previously published work. However, 
the subgap current is strongly suppressed upon the transition 
from the SIS to the SIN regime (right, B > 300 mT) and 
approaches the theoretically predicted limit. The inset shows 
an electron micrograph of the junction. 
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FIG. 2: (Color online) Ratios G sg ,min/GN versus normal 
conductance per geometric area Gn/A in SIN junctions; ex- 
perimental points are indicated with dots and triangles, the 
lines show theoretical predictions of Eq. (1). Inset: For a 



fixed Gat, we observe that G s 



is inversely proportional 



to A, thus exhibiting a linear dependency on T as suggested 
byEq. (1). 



for G sg set by Andreev reflection, we fabricated and 
investigated generic SIN Al/AlO^/Cu junctions. Here 
we achieved similarly small values of G sg , as shown in 
Fig. [2] Moreover, we observed a zero-bias conductance 
peak (ZBCP), a fingerprint of coherent diffusive Andreev 
transport [22-25:. This peak is suppressed by magnetic 
field and thus was not observed in the all-Al junctions. 
We further measured a dependence of G sg on the effec- 
tive electronic mean free path in the electrodes. As we 
will show, these observations indicate homogeneity of the 
tunnel barriers on the spatial scale of the mean free path. 

Our junctions were fabricated on oxidized Si substrates 
by the standard Dolan evaporation technique with in-situ 



thermal growth of aluminum oxide. Both metal leads are 
evaporated resistively at a base pressure of ~10 -6 mbar 
with a 0.5-1 nm/s evaporation rate. The junctions were 
e-beam patterned using a two or three layer resist sys- 
tem. The A1/A10 X /A1 junctions were ~100-400 nm wide 
with an overlap of ~400 nm [Fig. [ija)]. A variation of 
the oxygen dose (the product of pressure and time) from 
17 mbar • s to 8 bar • s yielded R^v values from 0.2 to 13 
kfi. The Al/AlCv/Cu junctions were ^80 nm wide with 
an overlap of ^220 nm. A variation of the dose from 
1.9 to 120 bar • s yielded Rat values from 1.5 to 200 kfi. 
We note that the growth conditions for the base elec- 
trode and barrier are the same in both types of devices. 
We therefore assume that the barrier qualities are simi- 
lar in both. The junctions were measured in a dilution 
refrigerator at temperatures below 100 mK. Each dc line 
was equipped with a two-stage RLC filter, a powder fil- 
ter, and with 2 m of thermocoax cable. The junctions 
were voltage biased and the current was read out either 
through a bias resistor or a transimpedance amplifier. 
The former method introduced excess noise for the high- 
Rjv junctions which is why the latter method was used 
for these junctions. 

To understand the experimental data, we applied a 
theory of the superconducting proximity effect based 
on quasiclassical Keldysh-Green function techniques [551 
|2"3"] . Both single- and two-particle processes contribute 
to the subgap conductance. The single-particle contribu- 
tion is given by the conventional tunnel model [4 . The 
two-particle Andreev conductance is given by 
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Here f N (E,V) = (l/2){tanh[(E+eV)/2k B T}-ta,Tih[(E- 
eV)/2k B T}} and E Th = T1D/L 2 is the Thouless energy 
with D the diffusion constant, / the mean free path, and L 
the length of the normal electrode. Eqn. ([I]) is valid in the 
experimentally relevant limit of phase-coherent diffusive 
transport over the distance L 3> I and assuming that 
Gjv is significantly smaller than the conductance of the 
normal lead. Furthermore, it is assumed that Etk "C A. 

Important features of the Andreev conductance 
which are relevant for interpretation of the experiment 
are i) the presence of an additional factor T indicating 
the two-particle origin of the Andreev transport, ii) the 
enhanced value at zero voltage, Ga(0) = GnT(3L/41) 
(ZBCP), and iii) the dependence on I. 

The ZBCP is indicative of the coherent, diffusive trans- 
port regime. It is explained by an electron-hole trans- 
mission resonance formed by the interplay between An- 
dreev reflection and scattering by impurities [23, 25]. The 
resonance is destroyed by an external magnetic field B, 
and the ZBCP disappears according to Ga(V = 0, B) = 
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FIG. 3: (Color online) (a) Subgap conductance vs voltage 
measured in SIN tunnel junctions with transparencies T = 
1 x 10~ 4 (blue, right y-axis) and T = 2 x 10~ 5 (red, left y- 
axis), solid lines. The best fit (dashed lines) is obtained by 
adding the single-particle conductance to Ga{V) of Eq. |T]) 
(dash-dotted lines), assuming an electronic temperature of 
100 mK, and finite values of the Dynes parameter (defined in 
the text), 7 = 4 x 10~ 5 A and 7 = 4 x 10~ 4 A, respectively, 
(b) Suppression of ZBCP by magnetic field B applied parallel 
to the substrate. Dots are experimental data. The line is 
theory. Inset shows G(V) curves for fields: B=0; 180 mT; 
280 mT (from top to bottom). 



G A {V = 0) tanh(&)/6, where b = ^/2\ L LeB/h, and A L is 
the London penetration length [22] . The ZBCP is usu- 
ally observed in high-transparency SIN junctions based 
on two-dimensional electron gases [23] , though it has also 
been found in metallic junctions with intermediate trans- 
parencies [XT] [2U] . Other transport characteristics of dif- 
fusive SIN junctions, e.g., shot noise, also exhibit a zero- 
bias anomaly [26) . 

Figure [2] shows the measured ratios G sg ^ m i n /GN for 
a wide range of barrier thicknesses in Al/AlO^/Al and 
Al/AlOz/Cu junctions measured in the SIN regime. Here 
G sg ,min indicates the minimum value of G sg approxi- 
mately at eV ps A/2. The conductance per unit area, 
Gn/A, where A is the junction area, was varied over 
more than 2 orders of magnitude. The measured ratios 



G 



I Gjv are proportional toG^/A as the theory (red 



sg,min 

and black lines) predicts, down to the lowest values mea- 
sured ~10 -5 . At a fixed Gn, we observe that G S g, m in 
is inversely proportional to A, indicating a linear depen- 
dence on T (inset of Fig. [2| . 

G sg vs voltage characteristics have been investigated 
in detail in all measured Al/AlO^/Cu junctions. Fig- 
ure [3] shows curves for two junctions, with parameters, 
G N = 50 nS,T = 2 X 1CT 5 (red solid), and G N = 
300 fiS,T = lx 10~ 4 (blue solid). In order to fit G N and 
the one- and two- particle contributions to G sg simulta- 
neously, one must make assumptions about the distribu- 
tion of r for the many microscopic conduction channels in 
the junction. We use a minimal model of a uniform T in 
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FIG. 4: (Color online) Scaling of the transport characteristics 
with thickness of the Cu electrode d: (a) specific resistance of 
the Cu electrode, (b) normalized zero- voltage resistance, and 
(c) normalized maximum value of resistance. Dots correspond 
to averages over several measured junctions, the error bars 
denote the standard deviation. Dashed lines show asymptotic 
behaviors extracted from Eq. |TJ. 



an active area which may be smaller than the geometric 
area. We find that the active area is about 13% of the ge- 
ometric area (cf. pj]). The other parameters used for the 
theoretical fitting are A = 210 /xeV, Dq u — 130 cm 2 /s, 
and Leu = 5 (tm. The best fit is achieved assuming an 
electronic temperature of 100 mK (the cryostat tempera- 
tures were 70 and 40 mK, respectively), and a tunneling 
density of states that is broadened due to spurious in- 
elastic processes. The broadening is parameterized by 
adding a small imaginary part 7 to the quasiparticle en- 
ergy, E+ij (the Dynes parameter [H][57]), having values 
of 7 = 4 x 10 _5 A and 7 = 4 x 10~ 4 A, respectively. Note 
that 7 affects the single-particle conductance but not Ga- 
We found similar values of temperature and 7 values of 
the same order for all junctions that were fit. 

We also measured the effect of an applied magnetic 
field on the ZBCP. The results [Fig. ^b)] are in a good 
quantitative agreement with the theory. The ZBCP dis- 
appears in magnetic fields larger than 280 mT, even while 
the Al electrode remains superconducting. 

To verify the dependence of the subgap conductance 
on the mean free path, we fabricated junctions with dif- 
ferent thicknesses d of the Cu layer. In bare Cu wires, 
resistivity as a function of d exhibited a crossover from a 
constant to a 1/d dependence at d ps 50 nm [Fig. Qa)], 
suggesting a crossover to surface dominated electron scat- 
tering. Figures Qb) and[4|c) show measurements of R sg 
vs d. Each point is the average of several junctions. At 
approximately the same thickness, d < 50 nm, the data 
exhibit a crossover from constant behavior to a variation 
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consistent with Eq. Q, namely R sg (V — 0) oc d and 

Rsg ,max C£ d. 

From the measured dependence of G sg on electrode 
thickness, we can draw important conclusions about the 
uniformity of our tunnel barriers. The employed theo- 
retical model assumes a homogeneous tunnel barrier on 
the scale of I. The fact that our observations agree with 
this model implies that any inhomogeneity of the tunnel 
barrier has a spatial scale larger than ~ 50 nm. If subgap 
transport were dominated by pinholes with a smaller size, 
the current would rapidly spread out in the electrode, so 
that electron scattering would not play a role, and trans- 
port would resemble that of a ballistic constriction. In 
that case, G sg would not show any dependence on d. We 
note that this 50 nm size scale matches that of the metal- 
lic grains in our base electrodes. This suggests a picture 
where the tunnel barrier is uniform on any given grain, 
but varies from grain to grain. The active tunneling area 
we extract then corresponds to a couple of grains. 

Having concluded that the tunnel barriers in our de- 
vices do not have pinholes, we must therefore conclude 
that the greatly enhanced subgap current when one-and- 
the-same device is measured in the SIS regime cannot 
be attributed to pinholes. It must be caused by other 
mechanisms, e.g., environmental resonances [28] . and it 
remains an open question. It is reasonable to extend this 
conclusion to other junctions fabricated using similar fab- 
rication techniques, which are in fact quite common. 

We note that 7 has implications for decoherence in su- 
perconducting qubits. In SIS junctions, it gives a residual 
conductance, G(V = 0) = Gnj/A, which should lead to 
relaxation. Considering, e.g., a transmon qubit, we ob- 
tain a relaxation time T\ = 7rA 2 /Huj 2 ^. Observed values 
of T\ ~ 6.5 /is [29] would imply 7/A ~ 10~ 4 , similar to 
our results, although this is not proof of causation. 

In conclusion, we have demonstrated a decrease of sub- 
gap current by 2 orders of magnitude in tunnel junctions 
as one of the superconducting electrodes is made normal. 
Good quantitative agreement with theory was observed 
in SIN junctions over a span of more than 2 orders of mag- 
nitude of the junction transparencies, T ~ 10~ 4 — 10~ 6 , 
with a minimum value of G sg /GN ~ 10 -5 . We observed 
all the generic features of coherent diffusive Andreev 
transport. Taken together, these observations strongly 
suggest that highly transparent, microscopic pinholes in 
the tunnel barrier are not the explanation for the ob- 
served excess subgap current in highly resistive SIS tun- 
nel junctions. 
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